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We characterize the nonhnear evolution of the baryon acoustic feature as traced by the dark matter and halos, 
using a combination of perturbation theory and N-body simulations. We confirm that the acoustic peak traced by 
the dark matter is both broadened and shifted as structure forms, and that this shift is well described by second- 
order perturbation theory. These shifts persist for dark matter halos, and are a simple function of halo bias, with 
the shift (mostly) increasing with increasing bias. Extending our perturbation theory results to halos with simple 
two parameter bias models (both in Lagrangian and Eulerian space) quantitatively explains the observed shifts. 
In particular, we demonstrate that there are additional terms that contribute to the shift that are absent for the 
matter. At z = for currently favored cosmologies, the matter shows shifts of ~ 0.5%, 6=1 halos shift the 
acoustic scale by ^ 0.2%, while 6 = 2 halos shift it by ~ 0.5%; these shifts decrease by the square of the 
growth factor D{z) at higher redshifts. These results are easily generalized to galaxies within the halo model, 
where we show that simple galaxy models show marginally larger shifts than the correspondingly biased halos, 
due to the contribution of satellites in high mass halos. While our focus here is on real space, our results make 
specific predictions for redshift space. For currently favored cosmological models, we find that the shifts for 
halos at z = increase by ^ 0.3%; at high z, they increase by ~ 0.5% . Our results demonstrate that these 
theoretical systematics are smaller than the statistical precision of upcoming surveys, even if one ignored the 
corrections discussed here. Simple modeling, along the lines discussed here, has the potential to reduce these 
systematics to below the levels of cosmic variance limited surveys. 

PACS numbers: 95.36.+x,98.80.-k 



I. INTRODUCTION 

It has been known for many years that the coupling of pho- 
tons and baryons in the early universe results in an almost 
harmonic series of oscillations in the matter power spectrum 
lE. QJstlwith a scale set by the sound horizon, s ~ 100 Mpc 
(see H Is! for a detailed description of the physics in modern 
cosmologies and for a comparison of Fourier and configu- 
ration space pictures). This feature can be used as a 'standard 
ruler' to measure the expansion rate of the Universe, and this 
baryon acoustic oscillation (BAO) method is an integral part 
of current and next-generation dark energy experiments. 

While the early Universe physics is linear and well un- 
derstood, the low redshift observations are complicated by 
the nonlinear evolution of matter and the non-trivial rela- 
tion between galaxies and dark matter. The nonlinear evo- 
lution leads to a damping of the oscillations on small scales 
il H, S H [M [OLand the generation of a small out-of- 
phase component | ^ % [13, [U IS fisll . The damping of the 
linear power spectrum (or equivalently the smoothing of the 
correlation function) reduces the contrast of the feature and 
thereby the precision with which the size of ruler may be mea- 
sured. The out-of -phase component corresponds to a shift in 
the acoustic scale which would bias the distance measure if it 
were not taken into account. 

In this paper we investigate the behavior of the acoustic sig- 
nal in biased tracers of the nonlinear mass field. We find that 
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biased tracers have different shifts than the matter, and discuss 
how these shifts can be modeled and significantly reduced. 

II. PRELIMINARIES 

One of the challenges of calibrating systematic effects in 
BAO is the very large scale of the acoustic feature, which 
requires huge volumes to be surveyed/simulated and reduces 
the effects of nonlinearities and astrophysical uncertainties. 
While the latter is what makes BAO an attractive standard 
ruler, the combination makes it challenging to measure sys- 
tematics with any statistical precision. To avoid these issues, 
we start with a toy cosmology (cCDM, see also lfl6ll ) which 
has iljv/ = 1, VIb = 0.4, h ~ 0.5, n = 1 and (78=1 • This cos- 
mology has an unrealistically high baryon fraction, a much 
smaller acoustic oscillation scale 50/i~^Mpc compared 
with the ^ 100 /i~^Mpc of the 'concordance' cosmology) and 
is more nonlinear at z = than our Universe is believed to 
be. This emphasizes the effects we are investigating while 
reducing the sampling eiTor, simplifying the numerical prob- 
lem and allowing us to obtain highly robust measures of small 
effects. In § |V] we extend the model constructed to concor- 
dance ACDM cosmologies, focusing on one with JIm = 0.25, 
fis/i^ = 0.0224, h = 0.72, n.s 0.97 and cts = 0.8 for defi- 
niteness. 

To model nonlinear structure formation and the formation 
of dark matter halos we used 10 independent simulations each 
of 1024^ particles in periodic, cubical boxes of side length 
2 /i~^Gpc. The simulations were started at z = 100 using 
the Zel'dovich approximation and evolved to z = with the 
TreePM UtIi code. The full phase-space data were dumped at 
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FIG. 1 : The mass functions at jz = for our cCDM (triangles) and 
ACDM (squares) cosmologies, along with two commonly used fit- 
ting functions due to Press & Schechter 03] (dotted, blue) and Sheth 
& Tormen IT^ (dashed, red). Since we are using the sum of the par- 
ticles in a FoF group for our definition of mass, we do not expect 
perfect agreement with either fitting function. The simulations for 
the ACDM cosmology are discussed in ijV] 
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FIG. 2: The nonlinear matter power spectrum at 2 = (red, up- 
per) and 2 = 1 (blue, middle), compared to the linear theory (black, 
lower) for our cCDM and ACDM models. The 2 = 1 power spectra 
have been scaled by 1/ to match the other power spectra on large 
scales and the ACDM spectra have been offset (vertically) for clar- 
ity. Note the strong damping of the oscillations and the large excess 
power on small scales in the evolved fields. 



a number of redshifts between 2 = 1 and 2 = 0, and groups 
were found using the friends-of-friends algorithm ifisll with a 
linking length of 0.168 times the mean interparticle spacing. 
We keep all groups down to 10 particles, or 2 x 10^^ H^^Mq, 
using the sum of the particle masses in the group as our halo 
mass definition for simplicity. These minimum masses corre- 
spond to peak heights running from i^~lat2 = 0toi^~2 
at 2 = 1. The mass function and nonlinear power spectrum 
for this model are shown in Figures [1] and |2] for reference. 
Throughout this paper we do not subtract shot-noise from the 
power spectra, but allow it to be a nuisance parameter in our 
fits (see below). More details on these simulations, including 
convergence tests, are in (l^ . 

Fitting the acoustic scale involves locating the position of 
the acoustic feature while allowing for variations in the broad- 
band shape due to e.g. galaxy biasing. We do so by fitting the 
observed power spectra with 

Pnt{k) = B{k)P^{k,a)+A{k) , (1) 

where A{k) and B{k) are smooth functions and a measures 
the acoustic scale relative to a "best-guess" cosmology. is 
a template for the biased, nonlinear acoustic feature. This def- 
inition of peak "shift" is over and above the shift in the point 
where ^'{r) = 0, or shifts in the extrema of the oscillations in 
the power spectrum. A good match between theory and ob- 
servations, including the correct background cosmology and 
hence distance-redshift relation, should give a = 1. 

Note that the precise partitioning into acoustic feature and 
broad band shape is dependent on the particular choices of A 
and B. Since constructing an accurate template for the acous- 



tic feature yields a good template for larger scales, we assume 
B{k) is a constant. A{k) is assumed to be a cubic spline spec- 
ified at 0.0, 0.1, • • •, 0.4 and derivatives specified at the end 
points. The shot noise component is simply absorbed into 
A{k). The above prescription yields seven nuisance param- 
eters; we do not vary this (or the particular prescription) in 
this paper, although we do get consistent results for different 
choices. 

The fits are done by minimization, fitting the 70 
power spectrum bins between k = 0.02 /i Mpc^^ and 
0.35/iMpc^^. We assume a diagonal covariance matrix 
where the errors are a smooth fit to the run to run variance 
of the 10 simulations. The errors on all derived quantities are 
determined the variance of 1,000 bootstrap resamplings. 

We now turn to the purpose of this paper - the determination 
of Pw, first for the matter and then for biased tracers. 



III. MATTER 

As is dramatically evident in Fig. |2] the sharp acoustic fea- 
ture at high redshift gets smeared by bulk flows and super- 
cluster formation as the Universe evolves IH 0, HI IH O . 
An estimate of the smearing is given by convolving the mat- 
ter correlation function with a Gaussian of width equal to the 
rms displacement of particles from their initial positions. To 
lowest order this is the rms Zel'dovich displacement 

^1 = A / dkPLik) (2) 
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FIG. 3: The shift in the acoustic scale, a — 1 vs. redshift for the mass 
(black triangles) and v — 1.9 halos (red squares) using Eq. ([3} as the 
template. Also shown are the best fit power-laws (dashed) and the 
expectations of perturbation theory [a — 1 oc D^(2)] (dotted). The 
points and curves for the halos are shifted by (Jz: = 0.05 for clarity. 



suggesting a template of the form 



Pw(fc, a) = exp 



PL(k/a) 



(3) 



where Pl is the linear theory power spectrum at the redshift 
of interest, and S is allowed to vary. The lowest order result is 
the same whether obtained using the peak-background split 
10], renormalized perturbation theory or resummed La- 
grangian perturbation theory |10:i]. While the Gaussian form 
of the smearing is a reasonable approximation to what we see 
in the simulations, we find 0(10%) deviations from Si in our 
cosmology. 

Figure[3]and TableUsummarize the shifts obtained using the 
above template. We find that a 7^ 1 at very high significance, 
with a — 1 ~ 4% at 2: = 0. We remind the reader that this is an 
extreme cosmology; the shifts for a concordance cosmology 
(as we discuss in f|V]i are approximately an order of magnitude 
smaller. We also observe that the shifts decrease with redshift 
as a — 1 ^ (adopting the convention D{z = 0) = 1), 
consistent with the field getting more linear at higher redshift 
and suggestive that the next order terms in perturbation theory 
(OiPl), see iH for a recent review) are responsible (see also 
iol llir who emphasized this point). 

To explore this possibility, we expand the density contrast 
in powers of linear density 5l 

6 = (5(1) + 5(2) + <5(3) + . . . (4) 

with S^^'' = It is straightforward to show that 
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TABLE I: Shifts, i.e. a — 1 in per cent, for the matter density au- 
topower spectrum and the cross-spectrum of the density with the lin- 
early extrapolated initial density field. The presence of a shift in the 
first column and not in the second demonstrates that the shifts arise 
from the higher order Pmn {m, n > 2) contributions to the power 
spectrum (see text for details). The last column demonstrates that 
this shift can be corrected by adding a P22 term to the template. 
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no. 4: The out-of-phase contribution predicted by perturbation the- 
ory well approximates the derivative of the acoustic signal. The 
points plot P22, while the line is a scaled version of (IPl /dink with 
the scaling given in the inset. The smooth components of both curves 
have been subtracted by fitting a cubic spline to the data. All curves 
are for z — 0. 



where the Fn contain dot products of the vectors qi and can be 
generated from recurrence relations lfl9ll20ll2lll . If the initial 
field is Gaussian the nonlinear power spectrum is then given 
by 

Pnl = {Pii + Pi3 ± Pi5 + • • ■} + {P22 + ■■■} (6) 

where Py (5(*)(5(j)) andPn Pl. In what follows, we re- 
fer to these groups of terms as the Pi„ and Pmn terms respec- 
tively. The Pin terms also arise if we take the cross-spectrum 
of the initial and evolved fields, while the Pmn terms only 
arise in the auto-spectrum. 

Ignoring combinatorial factors, the Pi„ terms are given by 

„("-l)/2 

Pln{k) PL{k) Yl [d^kPLiqk)] Fni- ■ ■) . (7) 

fc=l 

If we factor out the common Pl and focus on the lowest order 
correction, P13, we see that it involves a integral over a single 
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Pl and a relatively broad kernel which suppresses the oscil- 
lations. We expect these terms not to lead to significant out- 
of -phase contributions, though they can contribute corrections 
to the damping described above. Since we can isolate these 
terms by considering the cross-spectrum between the initial 
and evolved matter fields, we can test the above hypothesis. 
Table H] demonstrates that the shifts in this cross-spectrum are 
reduced by over an order of magnitude compared to the auto- 
spectrum and are consistent with zero given our statistical pre- 
cision. Note that the above argument is only true for the lowest 
order contribution, but higher order terms are suppressed by 
additional factors of 0((5^), and therefore drop off even more 
strongly with redshift. 

The lowest order P,„„ term is 

P22ik) = ^Qi{k) + ^Q2{k) + ^Q,{k) , (8) 

where the Qs are defined in Appendix |A] In constrast with 
Pi3, these terms (see Eq. lAlb involve integrals of products of 
Pl which can lead to out-of -phase terms. For example, F2 is 
peaked around « 52 ~ k/2. When Pl contains an oscil- 
latory piece, e.g. sin(fcr), P22 contains a piece schematically 
of the form sin^(fcr/2) ~ 1 + cos(fcr), which oscillates out- 
of -phase with Pl. Panel (a) of Fig. |6] explicitly shows this; 
in fact, this oscillatory part of P22 is very similar to scaled 
log-derivative of Pl (Fig.|4l see also i^). Taylor expanding a 
shifted power spectrum, 

Fi(fc/«)c.Pi(fc)-(a~l)^ + --- (9) 

we find good agreement between the predicted shift of Fig.|4] 
and the measured shift in Table H] It is important to note that 
we have subtracted smooth components for all of these com- 
parisons suggesting that even though perturbation theory does 
not accurately predict the broad-band shape (idi , it does cap- 
ture the evolution of the acoustic feature. 
The above suggests a modified template, 

Pw(fc, a) = cxp f — j Phik/a) 

+ expf P2i(kla). (10) 

Note that we have damped the oscillations in P22 as for P^ 
although we fix the damping scale to the first order calcula- 
tion. While this damping follows naturally from the heuristic 
picture described at the start of this section, it is also a con- 
sequence of the resummations in eg. Lagrangian perturbation 
theory. Table U shows that such a template corrects for the 
shifts observed in the matter. 

Looking ahead to biased tracers, we note that one cannot 
simultaneously fit for the amphtude of the P22 term and the 
shift, since these are highly degenerate (Eq. |9]l. Doing so re- 
sults in highly degraded constraints on the acoustic scale, as is 
evident in Fig. |5] We need to know the relative amplitude of 
P22 and Pii to correct the shift. For the mass the relative am- 
plitude is straightforwardly given by perturbation theory. Do 
the same terms come in for biased tracers and are we able to 
determine the relative amplitude of the two types of terms? 
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FIG. 5: The /S.y^ for fits of our z — mass power spectrum to 
the functional form of Eq. [5] (dashed red), to the form including P22 
(solid black) and marginalizing over the amplitude of the P22 term 
(dashed blue). There are 60 degrees of freedom in the fit, and in 
each case the best fit is a reasonable fit. Note that Eq. [3] gives a 
biased acoustic scale, including the P22 term eliminates the bias, and 
allowing the amplitude of the P22 term to float results in very weak 
constraints. 



IV. HALOS 



A. Shifts 



We investigate the acoustic signal of biased tracers in our 
simulations by computing the clustering of samples of dark 
matter halos chosen to lie in narrow mass ranges. Specifically 
we use the linear theory power spectrum to convert from halo 
mass to peak height, v = 5c/cf{M), and pick halos in the 
range 0.85 < v < 1.15, 1.15 < v < 1.45, • • •. For these we 
compute both the auto-power spectrum and the cross-power 
spectra with the linear and evolved dark matter density field. 
Given the low number density of most of our samples, we 
focus on the cross-power spectra in the analysis below. 

Motivated by the development in the previous section, we 
test if analogous results exist for halos; Table HI] and Fig. [5] 
summarize our findings. We find that (i) the halos exhibit 
non-zero shifts that are functions of halo type, (ii) the shifts 
scale approximately as D^, and (iii) the shifts are once again 
absent in the cross-spectrum with the linear density field. As 
with the matter, this argues that, within the language of pertur- 
bation theory, the shifts come from P,„„ terms and are domi- 
nated by second order corrections. The amplitude of the P„i„ 
terms relative to the Pn terms depends on the type of tracer, 
to which we now turn. 
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FIG. 6: The six combinations of Qn{k), including tlie exponential damping (at z = 0) that appear in the nonlinear halo power spectrum 
for cCDM. For each of these, we subtract a smooth component (by fitting a five point cubic spline), and compare to the "no-wiggle" power 
spectrum of iQ]. The dashed [red] line shows the same procedure applied to the linear power spectrum (divided by 5), while the dotted [blue] 
line is the P22 correction to the nonlinear matter power (upper left panel). In each panel we have scaled P22 by the multiplicative factor shown 
to better match each combination of Q„. 
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TABLE II: Shifts in the measured acoustic scale probed by the halo 
density field as a function of halo peak height and redshift. The first 
group of numbers are the measured shifts for the halo density field 
correlated with the matter field, fit with the template of Eq. [S] The 
second group are for the halo density correlated with the linear den- 
sity field, again fit with the template of Eq.[3] The final set of num- 
bers are analogous to the first, except fit using the template including 
the Pmn corrections. Note that the shifts present in the first group 
are significantly reduced in the other groups. 



B. Eulerian and Lagrangian Bias 

We can proceed to develop the perturbation theory of biased 
tracers in two ways: via Eulerian or Lagrangian perturbation 
theory. We begin with the former and follow ll22ll in defining 

SH^bfS+^-ls' + .... (11) 

Implicit here is that the halo density is defined in configuration 
space, and that the density fields have been smoothed on some 
scale R to allow us to truncate the expansion. We assume that 
we are working on scales ^ R, and will ignore subtleties 
that arise from the smoothing ll23l l24ll for now; we explicitly 
reinstate the smoothing scale in Sec. |IVC~2] 

The halo auto-power spectrum in Eulerian perturbation the- 
ory becomes (see also ll23ll24il25ll ) 

P„ = {bff{Pn+P22) + bfbi (^^Qs + Q9^ 

+ ^^Qi3 + --- (12) 

with terms like Pi„ included in the missing terms denoted 
• • •. The cross spectrum between two tracers / and // can be 
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obtained by the replacements 



bl 
6162 



^{1)^(11) 
"1 "2 



(13) 



As an example, for the cross-spectrum with the mass we ob- 
tain 



Ph,,n^bf (Pll +P22) 



6f 6f ^3^ 



(14) 



Figure|6]shows that the P22 and jQs + Qg terms contain out- 
of -phase oscillations with very similar shapes while the Q13 
term is essentially non-oscillatory. This suggests that biased 
tracers will exhibit different shifts than the matter, and the dif- 
ference will depend on the structure of the bias. It is also 
worth pointing out that the shift for a 5i = 1 tracer is not the 
same as for the mass - a fact also evident in Tablellllwhere the 
ly = 1 halos exhibit different shifts from the matter! 

The alternative description is within the Lagrangian pic- 
ture, which has recently been developed in ll26ll (see also Ap- 
pendix|B]for the basic definitions). Within this formahsm the 
halo auto-spectrum can be written 



Ph = 



exp 
bi 



P11+P22 

3 



Q5 + 2Q7 

{biy[Q9+Qii] 

26f&^Qi2+ ^(6^)'Qi3| 



Q8 + Q9 



(15) 



where again terms like Pi„ have been included in • • •. As 
before, expressions for cross-spectra follow from the mapping 
in Eq. [T3] taking care to expand the (1 + 6f )^ term before 
making the substitutions. Again, Figure|6]shows that the terms 
which arise look like scaled versions of P22, except for Q13 
which is non-oscillatory. 

The structure of these two sets of power spectra appears 
quite different, especially in the scaling of the different out- 
of -phase terms with 6„. However once broad-band power is 
removed, both of these cases can be effectively written as 



Ph 



exp 



[BiPl + 62P22] 



(16) 



where we have implicitly assumed that in both cases some of 
the higher order terms we have neglected above would sum 
to an exponential damping, as happens in some variants of 
both Eulerian and Lagrangian perturbation theory. For the au- 
topower spectra, using the empirically determed scaUngs in 
Figure |6l the Bs are related to the bias parameters by 
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(17) 
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FIG. 7: (Upper) The best fit values of B\ and Bi for the cross power 
spectra of halos and the evolved matter density; triangles, squares, 
circles, and crosses are the 2 = 0, 0.3, 0.7 and 1 data respectively. 
Note that increasin g B \ corresponds to increasing v. The solid lines 
are Sheth-Tormen 111 511 predictions, while the dashed lines are for 
Press-Schecter IT4II : the thick [red] lines are for Lagrangian theory, 
while the thin [blue] lines are for Eulerian theory. The dot-dashed 
lines are based on a quadratic fit to 62(61), again both for Eulerian 
and Lagrangian bias models. The theoretical scatter in these relations 
is ^ 20%. (Lower) As above but for the halo-halo auto-spectrum. 



or 



L\2 



Bl = (1 + fet 
B2 = ^ + \b{ + \bi 



b[b^2 



(18) 



Analogous expressions can be written for the cross-spectra by 
making the substitutions described above. 

If we fit our N-body data to Eq. ( fTSI ) we find that this form 
is a good description of the data and the different samples all 
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lie in a narrow band in the Bi — B2 plane, as shown in Figure 
|2]for both the auto- and cross-spectra. This suggests that the 
halos form a 1 -parameter family in terms of the nonlinear bias. 
This is fortunate, because the ratio of B2 to Bi is degenerate 
with the shift of the acoustic scale. 



and Lagrangian densities may be related by the Taylor series 
iH (5^ = (5^ + c{Se)^ with c ~ -0.805 assuming spherical 
collapse, we obtain (see also IS) 



6f = 1 + 6f 



(25) 



C. Explaining the Shifts 



and 



2fe[(l 



(26) 



The previous results simply imply that the amplitude of the 
shifts is a function of halo bias. However, the perturbative 
formulation of the previous section also relates the amplitude 
of the shift to the bias parameters of the halos. We test this 
relationship here, first using the peak-background split model 
(see e.g. |[l5l|23>|23l) and an empirically calibrated 61-62 re- 
lationship from simulations. 



1. Peak-Background Split 

The starting point for the peak-background split is the un- 
conditional multiplicity function 



■'^ ' 2p dM 



(19) 



The theory then makes predictions for both the mass-halo 
cross-spectrum and the halo-halo auto-spectrum, with the for- 
mer having significantly less shot noise. Fig.|2]compares these 
predictions to the observed relation between the Ss, where we 
have translated from the bias parameters using the prescrip- 
tion in the previous section. For the cross-spectrum the simple 
model describes the observed trend with Bi, with the differ- 
ences between the Press-Schecter and Sheth-Tormen predic- 
tions being small (and not distinguishable by the data). Both 
models overpredict B2 at high Bi. Interestingly, the predic- 
tions for the Eulerian and Lagrangian descriptions are virtu- 
ally indistinguishable, even though the structure of the ex- 
pressions for the power spectra are very different. For the 
auto-spectrum we have much more limited N-body data, but 
again the simple model does an adequate, if not perfect, job of 
describing what we see. 



which can be fit with 



vf{v) cx 1 + 



av 



1/2 



cxp 



-^1 (20. 



where a = 1, p = gives the Press-Schecter mass function 
lfl4ll . while a = 0.707, p = 0.3 yields the Sheth-Tormen mass 
function ifisll . Within the assumption of the peak-background 
split, the conditional multiplicity function is given by the sub- 
stitution, 



(21) 



where 5 is the background density and 5c — 1.686 is the crit- 
ical overdensity for collapse. The Lagrangian bias parame- 
ters then follow from Taylor expanding the (appropriately nor- 
malized) conditional multiplicity function as a function of 5, 
yielding = [v f {v)]"^ dJ" / d6''[v f {v)] or 



1 



2p 



1 + (ai/2)P 



and 



1 



2p{2av'^ + 2p - 1) 



1 + {av^Y 



(22) 



(23) 



The Eulerian bias parameters are then defined by the mapping 
of the halo density 5^ from Lagrangian to Eulerian space. 



(24) 



where the factor of ( 1 + 5^ ) comes from the mapping of the 
Lagrangian to Eulerian volumes. Assuming that the Eulerian 



2. An empirical 61 — 62 relation 

A second approach is to empirically calibrate 61 and 62 us- 
ing models based upon simulations. Recall that in any survey 
aiming to measure BAO there will be ample data on small 
scales with which to construct models of the tracers. 

The traditional approach to determining the Eulerian bias 
parameters is to compare various moments of "counts-in- 
cells" to the perturbative expressions. A concern with such 
an approach is the effective scale of the measurements, and 
the validity of the perturbative expressions. An alternative is 
to compare the auto-spectra of the halos and the mass with the 
cross-spectrum. The combination of the three spectra can be 
used to isolate 61 and 62, but this tends to be very noisy and 
proper shot-noise subtraction is an issue. We outline a differ- 
ent approach below that explicitly only uses large scales; we 
defer detailed comparisons with other methods to future work. 

We start by considering the configuration space statistic 

5(x) = (5i(xi)<5„(xi)) 

^ 6f(<5i(xi)<55(x2)) + f (5i(xi)5|(x2)) ,(27) 



where x = xi — X2 and 5s is the nonlinear matter density 
smoothed on a scale R such that the Eq.[TT]is valid. Working 
to second order in the density field, the second term above 
reduces to 62((5l(xi)(55(x2))^. If we work on large scales 
|xi — X2 1 ^ R, then we can approximate the smooth fields by 
the underlying density field. Fourier transforming, we obtain 



(28) 
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0.05 h Mpc ^ to Eq.|28]where 6f is determined from the low 
k limit of the halo-linear density cross correlation. 

Fig. [8] shows the measured bias parameters for our cCDM 
simulations, compared with the Eulerian and Lagrangian pre- 
dictions. Fig. |2] demonstrates that the observed relationship 
well describes the observed Bi — B2 correlation, and therefore 
shifts in the acoustic scale. It is important to emphasize these 
constraints on 62 are independent of the acoustic oscillations. 
In fact, the 62 constraint depends on the Q13 contribution to 
the power spectrum, which is irrelevant for BAO. 

In principle higher order measures or the observed cluster- 
ing on small(er) scales contain information about Bi for the 
sample of interest, and we have shown that improving our 
ability to model the higher order terms could bear dividends. 
Testing such methods is beyond our scope here, and we defer 
it to future work. 



D. Shifts, Corrected - A Template for Halos 



FIG. 8: Measurements of 62 vs. 61 for our cCDM simulations. The 
symbols are as in Fig.|7] The solid [red] and dashed [blue] lines are 
the peak-background split predictions for the Lagrangian and Eule- 
rian bias, assuming a Sheth-Tormen Jish mass function. The dotted 
line is a simple quadratic fit to the data. 



The fc ^ limit yields a direct measure of 61", 
5(fc^0) = 6fgi3(0) 



where 



Qi3(0) 



{2nf 



[PL{q)f = I 



(29) 



(30) 



Note that these expressions explicitly work at large scales, 
where the perturbative expansion is valid. While the above 
derivation was for Eulerian bias. Appendix |B] shows that the 
same limit yields 



5(fc^0) = fe^Qi3(0) 



(31) 



implying that 63 = fef". This is a different relationship that 
what we obtained within the peak-background split, reflecting 
the different assumptions made. 

The procedure for determining b2 from our simulations 
are straightforward, except for one subtlety if determining Sl 
from the initial particle data. Since the initial particle posi- 
tions are tied to a grid, there is an excess of power at the par- 
ticle Nyquist frequency, which can alias to lower frequencies 
when computing . To avoid this, we smooth the initial 
density field before squaring; this modifies Qi3(0) to 



n3,s 



(0) 



(2^ 



[PLiq)Wiq)r 



(32) 



where W{q) is the Fourier transform of the smoothing kernel 
(we adopt a Gaussian of comoving width 5 h^^Mpc). We then 
determine 62 by fitting the S{k) measurements below k < 



We are now in a position to construct a template for the 
BAO feature traced by halos. The key ingredient is a cali- 
brated Bi — B2 relationship; we assume the Sheth-Tormen 
form of Eulerian peaks-bias here. Assuming an estimate of 
the large-scale bias, this fixes B2/B1. We fit the observed 
power spectrum to 



P^{k,a) = bi 



cxp [ ) PL{k/a) 



k^Y.l\ 62 p , , 



(33) 



where hi and I] are fit parameters. Si is determined from lin- 
ear perturbation theory (Eq.|2]). 

We show that this procedure returns (almost) unbiased es- 
timates of the acoustic scale in Table |ll] These results come 
from the halo-mass cross-spectrum, which is significantly bet- 
ter determined than the halo auto-correlation function. The 
results from the auto-correlation function are consistent with 
the shift being corrected, but the errors are too large to al- 
low a meaningful constraint with the simulations we have. In 
principle one could obtain even more accurate constraints by 
modeling each of the perturbation theory terms separately, but 
this becomes more model dependent so we don't pursue this 
line here. 

These results allow us to identify sources of systematic er- 
rors in the BAO measurement and estimate their level. How- 
ever, before doing so, we first extend our results from cCDM 
to ACDM. 



V. IMPLICATIONS FOR ACDM 

In order to extend the results of the previous sections to 
a ACDM cosmology, it is useful to summarize the various 
components of our model and see how they generaUze to a 
different cosmology. 



9 




s 
c 



0.1 0.15 0.2 0.25 0.3 

k (h/Mpc) 

FIG. 9: As in Fig.g] except for ACDM. 



• The shifts in the matter were caused by the P22 piece of 
the power spectrum, which well approximated a scaled 
derivative of the linear power spectrum (after subtract- 
ing out the broad-band shape). The scaling was directly 
interpre table as the shift in the acoustic scale. Fig. |9] 
shows that this continues to hold for ACDM, with the 
z ~ shift predicted to be ~ 0.5%, in agreement with 
simulation results by lfl2ll and analytic arguments by 

lilli. 

• Extending these results to biased tracers generated ad- 
ditional shifts, sourced by terms whose oscillatory com- 
ponents resembled scaled versions of P22- This contin- 
ues to be true in ACDM with exactly the same scalings 
of P22', Fig.fTOlshows an example. This implies that the 
template of Eq. [33] as well as the relationship between 
the Bm and the bias parameters 6„ continues to hold for 
ACDM. 

• The final component was to demonstrate that simple 
models of halo bias indeed explained the resulting 
shifts. The first of these - a peaks-bias model - is mani- 
festly cosmology independent. The second attempted 
to empirically calibrate the bias parameters; Fig. [TT] 
shows a similar calibration for ACDM. Interestingly, 
we find a similar cosmology independence for the em- 
pirical calibrations, manifested in both the comparisons 
with cCDM and with the different redshifts for ACDM. 

• The above allows us to predict that for ACDM the shift 
is roughly an order of magnitude smaller than cCDM 
and is given by a — 1 ^ 0.5% x x B2/B1. 

We can now turn to the systematic error budget for BAO. 
The expression above gives the bias in the acoustic scale if 
one ignored the effects in this paper. Errors in the calibration 
of B2/B1 directly translate into an error in the acoustic scale; 



0.05 0.1 0.15 0.2 0.25 0.3 
k 

FIG. 10: An example of the out-of-phase components, as in Fig. |6] 
except for a ACDM cosmology. Note the reduced amplitude of the 
Pmn terms, indicating smaller shifts than for our toy cosmology. We 
also note that the factors that scale P22 to the other Qn combinations 
appear to be independent of the underlying cosmology. 



eg. a 20% error in B2/B1 (approximately how good our toy 
models are) would correspond to a bias of 0.1% x D^. 

Our analysis also demonstrates that the degree of shift is 
sensitive to the degree of nonlinearity, or amplitude of the 
power spectrum. The out-of-phase terms scale as one higher 
power of Pl than the linear terms. To the extent that the am- 
plitude is degenerate with a change in bias of the tracer, uncer- 
tainty in the amplitude leads to uncertainty in the shift. As an 
example, if the local slope of the B2 — Bi relation is (3 and we 
imagine holding the large-scale power fixed, a change in the 
amplitude SPl/Pl = £ will induce a change 5Bi/Bi = — e 
and (2 — /3)e in the P22 term in Eq. ( fTSI l. From Fig. [T] we see 
typical values of /3 ^ 1. Applying the same scaling between 
a and P22 as above this would lead to a shift in the acoustic 
scale of - 0.005 e for ACDM. Thus 10% knowledge of the 
amplitude of Pl would give < 0.1% uncertainty in a. 

Fig. fT2l summarizes the systematic error budget, and com- 
pares it to the observational error goals for Stage III and Stage 
IV experiments l3lll . None of these systematics are expected 
to be relevant for Stage III experiments, and are within a factor 
of a few of the requirements for Stage IV experiments. 



VI. GALAXIES 

The above results have focused on the case of halo samples 
of a single mass, but can be generalized to arbitrary combina- 
tions of halo samples. Of particular interest is the halo model 
for galaxies (for a review, see |32:]), that has been very suc- 
cessful in describing the large scale clustering of galaxies. 

The halo model assumes that all galaxies live in dark mat- 
ter halos and the probability of a particular galaxy occupying 



10 



1.5 



0.5 

WW 



-0.5 



'I I 

y= 1.09x=-2.77x+1.31 





bE, ^ + 1 



no. 11: As in Fig. [8] except for ACDM. These measurements were 
based on an additional set of simulations employing 1200'^ parti- 
cles in cubic boxes of side 1250/i^^Mpc. The symbols - triangles, 
squares, circles, stars and crosses correspond to z = 0, 0.3, 0.5, 0.7 
and 1.0 respectively. The dotted line shows the fit from Fig.[8]while 
the dot-dashed line is the quadratic fit in the inset. 



a given halo depends only on the halo mass. Under this as- 
sumption the large scale clustering of galaxies is determined 
by the clustering of the halos, weighted by the mean number 
of galaxies in each halo. If we denote the mean number of 
galaxies in a halo of mass M as N{AI), the galaxy power 
spectrum is 



gal 



En 



(34) 



where P^^-' is the halo cross-spectrum for masses i and j and 
the weights are determined by the halo mass function, nh{M), 
and N{M) via 



w,=nh{M,)N{M,)AM, 



(35) 



Substituting Eq. [16] for the halo power spectrum and noting 
that the exponential damping is (to a good approximation) in- 
dependent of halo mass, we find that Pgai retains the same 
structure. 



^'gai(fc) = exp - 



[Sl,galPL + S2,galP22] • (36) 



but with new coefficients 



'n.gal 



AA2 



dMini^{Mi)N(Mi) X 



dM2ni-,{M2)N{M2)BniAh,M2) , (37) 



FIG. 12: Estimates of the shift (if not corrected) as a function of halo 
bias and redshift for our ACDM cosmology. The width of the shaded 
regions denotes the error estimated from the difference between the 
Press-Schechter and Sheth-Tormen forms in the conversion between 
Bi and B2/B1 and demonstrates approximately how errors in the 
theory propagate into a residual shift. We highlight two example pop- 
ulations of halos: a h{z = 0) = 1 sample [blue crosses] to represent 
emission line galaxies, and a h{z = 0) = 1.5 sample [red circles] 
to represent an elliptical sample. In both cases, the clustering of the 
sample is assumed to be constant with redshift. The errorbars corre- 
spond to a 10% measurement of the bias. Also shown are nominal 
distance accuracies for Stage III (currently underway) and Stage IV 
(future) experiments. 



where the normalizing factor 



N = dMnh{M)N{M) 



(38) 



is just the mean number of galaxies. Then, assuming 
S„(A/i, A/2) factorizes (or is a sum of factorizable pieces), 
we can simply replace 6„ in our earlier expressions with 



&n.gai = / dMnn{M)N{M)b,,{M) 



(39) 



Note that, as expected, Si,gai simplifies to fe^^j where the 
galaxy bias 6gai is just the weighted sum of the halo 61 's. The 
expressions for cross-spectra of galaxy samples follow a sim- 
ilar pattern to the halo cross spectra discussed earlier (Eq.[T3]). 
For a concrete example, we assume a TV (A/) of the form, 



N{M) = e (Af - A//„iin) [1 + A//Afi 



(40) 



We consider two cases : Mi ~ 00 ox & "threshold" sample, 
and Ml ~ lOA/min or a "satellite" sample; Fig. [T3]plots the 
shifts as we vary Mmin from IO^^/i^^A/q upwards. As one 
might expected, the increased weighting towards higher halo 
masses increases the shifts, although they remain smaller than 
the required systematics levels for near future surveys. The 
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FIG. 13: The shifts for three example galaxy samples at jz = for 
ACDM, as a function of the galaxy bias. The dotted line shows the 
shifts for halos, while the dashed [red] and solid [blue] lines are for 
threshold and satellite samples respectively (see the text for more 
details). We assume the Eulerian Sheth-Tormen peaks-bias model 
for definiteness. 



above also suggests that one could control systematics by re- 
weighting galaxies, the details of which will, of course, be 
population dependent. 



VII. REDSHIFT SPACE 

Until now, everything we have done has been in real space, 
ignoring the redshift space distortions that arise from pecu- 
liar velocities. The Lagrangian perturbation theory formalism 
allows us to predict the effect of large-scale redshift space dis- 
tortions in a straightforward manner, although the comparison 
with simulations is made more difficult and additional model- 
ing of the observations is required. We defer a detailed com- 
parison with simulations to future work, but comment on the 
trends here. 

The power spectrum now becomes anisotropic, P{k, fj), 
with ^ the cosine of the angle between the line-of-sight and 
k and our previous results correspond to /i = 0. Including 
the full /i dependence results in different combinations of Qn 
entering the expression. Some of these are in phase with Pl, 
while some are out of phase. This means the shifts in redshift 
space will be different than those in real space, as seen in sim- 
ulations (e.g. lfT2ll ). The in-phase term (£^12, see Eq. A73 in 
ll26ll ) is much smaller than P^, and we ignore it for now. 

If we concentrate on the isotropic piece of the power spec- 
trum JsH, we again find that the remaining combinations of 
Qn which enter are proportional to P22 and the constants 
of proportionality appear to cosmology independent (though 
they do depend additionally on / = dh\D/dhia ~ il"-^). 



There are no new degrees of freedom introduced theoretically, 
so in principle the redshift space shifts are determined from 
the same modeling as the real space shifts discussed previ- 
ously. 

For ACDM at z = the predicted shift in the matter 
grows from 0.5% to 0.75%, consistent with the shifts seen 
in lfl2ll . For biased tracers, the effect is to increase the shift by 
an /-dependent (but roughly bias-independent) constant. For 
/ = 1 (corresponding to high z), this constant is ~ 0.5% D^, 
while for / ~ 0.5 (corresponding to z ~ 0), it is ~ 0.3%. If 
uncoiTected, these shifts could be relevant for future experi- 
ments; we leave detailed calibrations to future work. 



VIII. DISCUSSION 

The propagation of acoustic waves in the early universe 
provides a robust means for determining the expansion his- 
tory of the Universe, and contraining cosmology. The stan- 
dard ruler is calibrated in the linear regime by observations 
of the CMB, but observed today with biased, nonlinear trac- 
ers. Since the acoustic scale is so large, the effects of bias and 
nonlinearity on the acoustic scale is small, but future exper- 
iments may have enough sensitivity that it needs to be taken 
into account. We have begun this program here. 

Using a set of N-body simulations of an extreme cosmol- 
ogy, in which the acoustic scale is relatively small and the 
nonlinearity quite pronounced, we have shown that shifts in 
the scale grow quadratically with the amplitude of the linear 
theory power spectrum for both the mass and for dark mat- 
ter halos. Motivated by this, and guided by arguments from 
Eulerian and Lagrangian perturbation theory, we found the 
dominant second-order contribution to the peak shift. This 
contribution, P22, quite well approximates the derivative of 
the acoustic signal, explaining why it leads to a peak shift and 
allowing us to estimate how the amplitude of the P22 term can 
be translated into a shift in the fitted acoustic scale. 

For dark matter halos the contribution depends on two bias 
parameters, bi and 62, allowing in principle arbitrary shifts of 
the acoustic scale. We showed however that dark matter ha- 
los, which will be the hosts of any galaxies we observe, obey 
a relation between bi and 62 which is relatively well predicted 
by the peak-background split. Once these two terms are re- 
lated the acoustic scale for nonlinear, biased tracers can be 
accurately determined, allowing high fidelity measurements 
of distances using baryon acoustic oscillations. 

We have described how redshift space distortions affect the 
scale shifts within the context of Lagrangian perturbation the- 
ory, where they increase the shift by tens of percent at low z 
and about a factor of 2 at high z. 

In this paper we have concentrated on the effects of non- 
linear gravitational evolution and halo biasing, showing that 
these effects could be understood at the 0.1% level. At this 
level of precision the inclusion of additional physics, such as 
differential evolution of the baryonic and dark matter compo- 
nents or the details of galaxy formation, may enter and in- 
vestigations in this direction should be undertaken. We have 
also assumed that the linear theory template is perfectly under- 
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stood, which also deserves further study. Such investigations 
must be part of any future high-precision BAO experiment. 

The simulations presented in this paper were carried out us- 
ing computing resources of the National Energy Research Sci- 
entific Computing Center and the Laboratory Research Com- 
puting project at Lawrence Berkeley National Laboratory. NP 
is supported by NASA HST-HF-0 1200.01 and LBNL. MW 
is supported by NASA and the DoE. This research was ad- 
ditionally supported by the Laboratory Directed Research 
and Development program at Lawrence Berkeley National 
Laboratory, and by the Director, Office of Science, of the 
U.S. Department of Energy under Contract No. DE-AC02- 
05CH11231. 



APPENDIX A: EXPRESSIONS FOR Pmn FOR m = 2, n = 2 

The expressions for P,„„ for a biased tracer (in both Eule- 
rian and Lagrangian perturbation theory) can be simply writ- 
ten by defining ||26i| 

= ^ drPL{kr) / dx PL{k^)Qn{r,x) , 

(Al) 

wherey(r,a;) = l + r^ — 2ra;, Pl is the linear power spectrum 
and the Qn are given by 



r2(l - x'^ f 

y 

— rxY 
rx(\ — x'^) 
X [1 — rxj 



(1 — x'^)rx{l — rx) 



\~x' 



y 

rx{\ — rx) 

y 



(1 -3ra-)(l 
V 



, Qio = 1 - .T , 



Qii = X , Qi2 = rx, Qi3 = r 



APPENDIX B: COMPUTING (^E^obj) IN RESUMMED LPT 

We present expressions for evaluating ([i5£](5obj) in re- 
summed Lagrangian perturbation theory, where 5obj is the 
density field of biased tracers. These expressions follow ll26ll 
and we refer the reader there for detailed calculations. 

The density field for a biased tracer can be defined by the 
displacement field \I'(q) and a function of the smoothed initial 
density field in Lagrangian space, F[6l{(1)], ^s 



Sohj (x) 



yA^[5i(q)]5gHx-q-*) , (Bl) 



where x and q are the Eulerian and Lagrangian positions and 
jj^' is the 3D Dirac S function. We implicitly assume that 
the argument to F has been smoothed on some scale much 
smaller than the large scales of relevance here, allowing us 
to ignore the smoothing here (see ll26ll for a detailed jus- 
tification). We cross correlate this with a field defined by 
cxp(iAi5L); (5£ is then simply obtained by taking the n-th 
derivative with respect to A and setting A to zero. 

The cross-power spectrum of the two fields is then given by 
(compare to Eq. 9 of 1.26!]) 



H(fc) 



d^qe 



— ikq 



ZTT 



i(Ai5i,(qi)+A2<5i,(q2))+ik*(q2) 



(B2) 



where q = qi — q2 and F is the Fourier transform of F. The 
correlators of interest are then given by 



([^El'^obj) = - 



1 d'^H 



dXl 



(B3) 



Ai=0 



The algebra now follows through as in 112611 using the cumulant 
expansion theorem, and collecting aU zero-lag correlators to 
yield (compare to Eq. 24 in ll26ll ). 



TLi^k) ~ exp 



m— 1 ^ ' 



exp 



d^qe-'^"^ 



dX2 
2tt 



FiX2) X 



-AiA2e(iqi)+ J2 E 



°° °° jTll+Tl2+ni2 



\ ^2 ^0m2 1^'^) 



rii +n2 > 1 ?7i2 > 1 

I 



(B4) 



where with (• • ■)c denoting the connected moments. 

Given Eq. IB4I it is straightforward (if tedious) to compute 
^{\^\) — {^l{<Ii)Sl{(12)): = C(0) , (B5) expressions for the correlators in Eq. lB3l Of particular interest 

to us here is the fc limit of ([<5£](5obj) which, as in Eulerian 
and perturbation theory, is given by 



pnin2 
-°0m2 



:[^L(qi)]"M<5L(qi)]"MkVI>(q2)]™^), , (B6) 



S{k ^ 0) = 6^Qi3(0) , 



(B7) 
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with 62 defined by 



d\ 
— f 
27r 



1 



dSe-'V^^^^ (B8) 
While the full formalism is required in general, we note that 



this fc limit can be obtained more simply by dropping 
the k • ^ term in Eq. IB2I leaving only Gaussian fields in the 
exponent. These fields have only a second connected moment, 
thus only the AiA2^(|q|) term survives in the last exponential 
of Eq. lB4l and expanding this exponential gives Eq.| 
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